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ABSTRACT 

Massive relic galaxies formed the bulk of their stellar component before z ~ 2 and 
have remained unaltered since then. Therefore, they represent a unique opportunity 
to study in great detail the frozen stellar population properties of those galaxies that 
populated the primitive Universe. We have combined optical to near-infrared line- 
strength indices in order to infer, out to 1.5 Re, the IMF of the nearby relic massive 
galaxy NGC 1277. The IMF of this galaxy is bottom-heavy at all radii, with the 
fraction of low-mass stars being at least a factor of two larger than that found in the 
Milky Way. The excess of low-mass stars is present throughout the galaxy, while the 
velocity dispersion profile shows a strong decrease with radius. This behaviour suggests 
that local velocity dispersion is not the only driver of the observed IMF variations seen 
among nearby early-type galaxies. In addition, the excess of low-mass stars shown in 
NGC 1277 could reflect the effect on the IMF of dramatically different and intense 
star formation processes at z 2, compared to the less extreme conditions observed 
in the local Universe. 
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1 INTRODUCTION 

A combination of major and minor dry mergers are the 
favoured explanation for the evolution of the scaling re¬ 
lations of massive galaxi e s since z ^ 2 (e.g. lNaab_^£_ab 
20091 : iTruiillo et al.l 201 ll : iFerreras et al.l lAl4 I Ruiz et~ 
20141 . Being a stochastic process, a natural prediction of 
this evolutionary channel is the presence, in the local vol¬ 
ume, of a number of massive galaxies which have experi¬ 
enced fe w if any significant m erger events since they were 
formed (lOuilis fc Truiillcll2013ll . These putative relic galax¬ 
ies present a unique opportunity to study in detail the prop¬ 
erties of massive, z ~ 2 galaxy analogues in the nearby Uni¬ 
verse. 

The Perseus cluster galaxy NGC 1277 fulfills all the 
characteristics to b e considered a relic massive galaxy 
dTruiillo et al.l[^14h . The galaxy exhibits an unusually high 
central velocity dispersion (above 400 km s“^) and is phys¬ 
ically com£a£tj_Mth_au effective circularized radius of 1.2 
kpc jTruiillo et al.ll2014h . The bulk of its stars are old (age 
> 10 Gyr) and a-enhanced ([Mg/Fe] ~ 0.3 dex), which sug¬ 


gests a short (r ~ 100 Myr) and intense (SFR > 10^ M 0 /yr) 
star formation burst at z ^ 2 , follo wed by subsequent quies¬ 
cent evolution (ITruiillo et al.l[2014l l. In addition, NGC 1277 
hosts the most massive ( compared to the total galax y mass) 
black hole foun d to date dvan den Bosch et al.ll2012l . but see 
lEmsellernll2013fl . 

ITruiillo et al.l (l2014h studied the stellar population 
properties of NGC 1277 assuming that its stellar init ial mass 
function (IMF) follow s that seen in the Milky Way (iKroupal 
I 2 OO 2 I : IChabrierll2003l l. Although this approach has been ex- 
tensi v ely used in the literature, man y studies ([Cenarro et al.l 
200^;_^n_^okkum&Conjo^[j2010klCa£2el]ari_et_^lJj20lF 


Ferreras et al. 20131 : 1^ Barbara et al.l '2013l : ^iniello et al.l 

2014h have shown that the shape and the normaliza- 


tion of the IMF cor r elate with galaxy mass. Moreover, 
iMarti'n-Navarro et al.l (l2015li have found that the IMF 
slope of some representative massive galaxies depends on 
galactocentric distance. Since the inferred properties of a 
stellar population stron gly depend on the adopted IMF 
llFerre-Mateu et ani2013ll . it is important to investigate the 
stellar populations of massive systems such as NGC 1277 
while relaxing the assumption of a universal IMF. 

The importance of measuring the behaviour of the IMF 
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as a function of radius in NGC 1277 is two-fold. On the 
one hand, it provides a route to understanding the physical 
conditions of massive galaxies when the majority of their 
stars were formed, with no influence from later accretion 
events. On the other hand, NGC 1277 shows an atypical 
combination of kinematical and stellar populatio n properties 
dvan den Bosch et al ] |2012l: iTruiilio et al. I I 2 OI 4 I I. which can 
be used to investigate the main driver behind these IMF 
variations. 

By studying the radial variation of IMF-sensitive spec¬ 
tral indices in NGC 1277, we have been able to derive its 
IMF radial profile, which remains bottom-heavy up to ~ 
1.5Re. We have simultaneously inferred, in a self-consistent 
manner (i.e., by adopting the best-fitting IMF value at each 
radius), its age, metallicity and [Mg/Fe] gradient s. Here we 
find that the claims made bv lTruiillo et al.l (l2014l l regarding 
NGC 1277 as being a relic object are also supported under 
the assumption of a steep IMF. 

This paper is structured as follows: in Section 2 we 
present our data; Section 3 contains the analysis and the 
main results, which are finally discussed in Section 4. We 
adopt a standard cosmology: Ho= 70 km s“^ Mpc“^, flm = 
0.3, and Ha=0.7. 


2 DATA AND DATA REDUCTION 


Our data consist of two sets of long-slit spectroscopic ob¬ 
servations of NGC 1277 carried out at the 4.2m William 
Hershel Telescope (WHT) and at the 10.4m Gran Telesco- 
pio Canarias (GTC), located at the Spanish Observatorio 
del Roque de los Muchachos on La Palma. The WHT spec¬ 
tra were acquired with the Intermediate dispersion Spectro¬ 
graph and Imaging System (ISIS), using the R300B grat¬ 
ing and a 1 arcsec slit, placed along the major axis of the 
galaxy fsee iTruiillo et aril2014l . for details). The total on- 
source integration time for these observations was 3 hours, 
with 1 arcsec seeing. The spectra cover the wavelength range 
from ~ 3700 to 6000 A, with 3.4 A spectral resolution 
(FWHM). In a second run of observations, we again targeted 
the NGC 1277 major axis with the GTC Optical System for 
Imaging and low-Intermediate-Resolution Integrated Spec¬ 
troscopy (OSIRIS) in order to increase the wavelength cover¬ 
age up to 10000 A, allowing us to inclu de near-infrared IMF - 
sensitive features such a s Na 8190 (ISch iavqn et al.l l200d l 
and the Calcium triplet dCenarro et al. 200 ih . These data 
have a (wavelength-dependent) resolution ranging from 6 A 
at A ~ 5000 A to 8 A at A 10000 A. The total exposure 
time was 1.5 hours on source with a typical 0.8 arcsec seeing. 
In both datasets, the pixel size was 0.25 arcsec. In this work, 
all measurements (i.e. line-strengths indices) blueward and 
redward of 5700 A come from WHT and GTC data, respec¬ 
tively. 

Dat a reduction wa s performed using the REDUCEME 
package (ICardiellll999l L which allows for the careful treat¬ 
ment of error propagation. The reduction process included 
the standard bias subtraction, flat-fielding, cosmic ray clean¬ 
ing, distortion (C and S) correction, wavelength calibration, 
sky subtraction and flux calibration. 

Note that, for our purposes, sky removal is a critical 
step of the reduction, as both airglow emission and telluric 
absorption have a prominent role in the raw spectra above 


A ~ 7000 A. To ensure an accurate sky subtraction, an 
on-and-off observational strategy was followed, where the 
galaxy was offset in the CCD plane during the observations. 
This allowed for a clean sky emission removal by subtract¬ 
ing consecutive exposures. The second crucial aspect was 
that of correcting the spectra for telluric absorption fea¬ 
tures. In order to investigate the impact of such correction 
on the derived stellar population parameters, we adopted 
two independent approaches. The spectrophotometric stan¬ 
dard star (Hiltner 600, Bl) was targeted at the end of the 
observations. Using the IRAF task telluric, the derived at¬ 
mospheric transmission spectrum was scaled to correct the 
NGC 1277 data. In addition, w e made use of the ESO tool 
Molecfit dKausch et al.l[2014| j. The latter does not require 
any extra calibration data (i.e. a spectrophotometric stan¬ 
dard star) as it constructs a synthetic atmospheric absorp¬ 
tion model for a given spectrum, by fitting spectral regions 
dominated by telluric absorption in the spectrum itself. 
To this effe ct, molecfit relies o n the radiative transfer code 
LBLRTM dClough et al.l[2005h . allowing for a simultaneous 
fit of both the telluric model and the line spread function 
to the data. The differences in the inferred stellar popula¬ 
tion properties among the two approaches are described in 
Appendix A. 

Since the analysis of stellar populations requires a 
higher signal-to-noise (SN) than kinematics, we performed a 
first radial binning by imposing a SN threshold of 20. The ra¬ 
dial velocity and velocity dispersion gradients, derived using 
the p ublic available software pPXF lICaPDellari fc EmsellemI 
l2004h . are shown in Fig. [T] Our meas urements are in good 
agreem ent with those published by Ivan den Bosch et al.l 
(I2OI2II . showing that NGC 1277 is pressure-supported in 
the central regions (ctc ~ 400 km s“^), whereas ordered 
motions dominate beyond ~ 0.5Re. We used the inferred 
kinematics (radial velocity and velocity dispersion) to cor¬ 
rect every row (spectrum) in the image to the rest-frame 
via linear interpolation of our measurements. Once all the 
spectra were at rest-frame, we calculated a representative 
spectrum at each radius of the galaxy with a minimum SN 
per A (at AA 6000,6200 A) of 80 in the outskirts, achieving 
over 700 in the central bin. 

Throughout this work, we assume a circularized effec¬ 
tive radius of Re = 1.2 kpc for NGC 1277, with a 0.344 kpc 
arcsec“^ scale (ITruiillo et al.lf^014^ . Fig. [5] shows three rep¬ 
resentative spectra at 0, 0.5, and 1 Re, respectively. The 
change in spectral resolution from the WHT to the OSIRIS 
data appears at A = 5700 A. Note also how the intrinsic (as¬ 
sociated with the velocity dispersion) resolution varies from 
the central (in red, a ~ 400 km s“^) to the outermost bin 
(in orange, a ~ 200 km s“^). A zoom around the H/3 feature 
is shown as an inset, suggesting no significant contamination 
of the line from nebular emission (see § |3I ■ 


^ In particular, we run Molecfit by fitting the spectral re¬ 
gions AA 8130,8300 A; AA 7585,7698 A; AA 7235,7326 A and 
A A 6861,6924 A, where strong H 2 O and O 2 telluric features are 
found. 
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Figure 2. Averaged spectra at 0.0 (red), 0.5 (blue) and 1.0 (yellow) Re, arbitrarily scaled in the vertical direction. Note the different 
spectral resolution between the GTC-OSIRIS (A > 5700 A) and the WHT (A < 5700 A) data. Note also the large variation in the velocity 
dispersion, reffected in a gradual smoothing when moving from the outer l.ORe spectrum to the central radial bin. The inset panel shows 
a zoom in of the H/3 line. 
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Figure 1. Velocity dispersion (solid line, left vertical axis) and 
radial velocity (dashed line, right vertical axis) profiles as a func¬ 
tion of the major-axis radial distance, normalized to the (circu¬ 
larized) effective radius Re. Dashed vertical line indicates the see¬ 
ing radius, measured by its average half width at half maximum 
(HWHM) value (0.4 arcsec). 


3 ANALYSIS AND RESULTS 

We based our stellar populati on analysis on the extende d 
version of the MILES models llVazdekis et al.ll20ld . |2012|) . 
The MILES models cover a range from —2.32 dex to +0.22 
dex in total metallicity, and from 0.06 Gyr to 17 Gyr in age. 
We assu med a bimodal IMF, p arametrised as a two-segment 
function llVazdekis et aL|[l996l ). where the star number per 
logarithmic mass bin is given by: 


We allowed for a variation in the IMF slope, Eb, from 0.8 
(bottom-light) to 3.3 (bottom-heavy). Under this parametri- 
sation, the standard Kroupa-like IMF is well represented by 
a slope of Fb = 1.3. Note that Fb varies the slope of the high- 
mass end of the (bimodal) IMF, in contrast to other studies 
where the high-mass end of t he IMF is kept hxed, while th e 
low-mass end is varied (e.g. iGonrov fc van Dokkurnll2012l ). 
Since the (bimodal) IMF is normalized to have a total mass 
of 1 Mq, changing Fb does actually vary the mass fraction 
of low- to high-mass stars. Hence, for the present study, 
Fb should not be regarded as the high-mass end slope of 
the IMF, but rather as a proxy for the fraction of low-mass 
stars. 

To determine the stellar pop ulation parameters we fo l- 
lowed the same approach as in iLa Barbera et ahl (I2ni3h . 
Age, metallicity, and IMF slope were obtained by minimis¬ 
ing the following equation: 
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X^(rb,age, [Z/H]) = ^ 

i 


{EWi 


Ac,,i) — EWM,i 

O'EWi 


( 1 ) 


where the sum expands over the ith observed and pre¬ 
dicted line-strength indices included in the analysis. EWi 
are the measured line-strengths, EWM,i the predictions from 
MILES models, and Aa,i is a correction applied to each in¬ 
dex to account for deviations from a non-solar abundance 
scale, and the effect of “residual” abu ndance ratios on line 
strengths Isee lLa Barbera et al.l[201^ and Appendix A, for 
details). To derive [Mg/Fe], for each radial bin we estimated 
separately the Magnesium and the Iron metallicities (ZMg 
and Zpe, respectively) through index-index fitting of the 

vs. Mgb and (Fe) line-strengths, respectively. We then 
subtracted one metallicity from the other to calculate the 
[ZMg/ZFe] proxy, whi ch tightly correlates w ith the “true” 
[Mg/Fe] (see Fig.6 in iLa I^rbera et ah 1 1201311 . In Fig. El we 
show the measured radial profiles of all indices considered 
in the fitting process, along with the best-fitting SSP so¬ 
lution at each radius. Notice that the line has been 
corrected for nebular emission contamination, as detailed in 
Appendix A. 

The derived age, metallicity and [Mg/Fe] radial pro¬ 
files are shown in Fig. [4] NGC 1277 exhibits a flat age pro¬ 
file, always above 10 Gyrs. The metallicity, on the contrary, 
decreases with increasing radius. Notice that the metallic¬ 
ity in the central bin (-1-0.41 dex) exceeds the maximum 
value of MILES models (-1-0.22). Hence, for this specific 
bin, we extrapolated t he models (in the same way as in 
iLa Barbera et aLll2013h to match the observed indices (in 
particular, the total metallicity indicator [MgFe]'). Although 
in general there is a good agreement between the observed 
and best-fitting line strengths in Fig. El some discrepan¬ 
cies exist. In particular, the mismatch between observed and 
model indices in the most central bins (see, e.g., panels (f) 
and (i)), can be at least partly understood as a limitation 
of the stellar population models to describe very metal-rich 
(> 0.2 dex) populations and/or uncertainties in the age es¬ 
timate. However, none of these issues significantly affect the 
results of our work, as discussed in Appendix A. Fig. E] also 
shows that the [Mg/Fe] gradient of NGC 1277 is flat, with 
[Mg/Fe]'~ 0.3 dex at all radii. This fiat [Mg/Fe] radial pro¬ 
file can be also directly inferred from the fact that both the 
Mgb and (Fe) spectral indices (see panels (a) and (b) of 
Fig. El show a similar amount of radial variation 0. Note 
that the large error bars on [Mg/Fe] in the central regions 
reflect the uncertainties on the [Fe/H] and [Mg/H] determi¬ 
nation/extrapolation used to derive the [Mg/Fe] proxy. The 
gradients shown in Fig . E] are in agreement with those pub¬ 
lished bv lTruiillo et al.l (l2014l . Fig. 4 in their paper). This is 
not a trivial result, since here we have allowed for possible 
radial variations of the IMF. 

The IMF radial profile of NGC 1277 is shown in Fig. El 
where the shaded region accounts for the scatter among 
AFb values from different tests, where the fitting process 


^ Notice, in fact, that in MILES models, both Mgb and (Fe) have 
a similar sensitivity to total metallicity, hence a similar radial 
gradient for the two indices implies a fiat trend of the [Zj^g/Zpa] 
proxy. 
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Figure 3. Observed gradients of the age-, metallicity- and IMF- 
sensitive indices considered in the fitting process (blue symbols). 
Dark blue represents the mean radial profile, whereas in light blue 
we show the radial gradient at each side of the galaxy. The Icr 
error bars are also shown for the mean profile. The best-fitting 
solution is overplo tted in orange. The def inition of all spectral 
indices are given in lLa Barbera et al.l ll201^b while aTiO has bee n 
computed according to the definition of rSpiniello et al. I lioll). 
Although not individually included in the fitting process, panels 
(a) and (b) show the Mgb and (Fe) line-strength profiles. Note 
that, according to the MILES models, both indices have quite 
similar sensitivities to a variation in total metallicity. Since they 
also exhibit the same radial variation O.SA), a fiat [Mg/Fe] 
gradient is expected, consistent with Fig. [4] In all panels, dashed 
vertical lines mark the seeing radius. 


has been repeated with different sets of spectral indices and 
different methods to remove telluric lines as detailed in Ap¬ 
pendix A. Regarding the discrepancies between the observed 
and best-fitting line strengths (see Fig. o, note that the 
Ti02 and NaI8190 lines can be potentially affected by er¬ 
rors in the flux calibration (the former) and telluric residu¬ 
als (both). Although the central Ti02 value might suggest a 
slightly steeper IMF proHle, our best-htting solution points 
to a mild radial variation of IMF slope (Fb ~ —0.5), i.e. a 
bottom-heavy IMF at all radii. Notice that in the present 
analysis we have not included the CaHl^ and CaH2 IMF- 

® Although not used in our analysis, the radial variation of the 
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Figure 4. Best-fitting age (top), metallicity (middle) and 
[Mg/Fe] (bottom) gradients. Both age and [Mg/Fe] show a fairly 
flat radial behaviour, whereas metallicity decreases by 0.4 dex 
from the centre outwards to 1.5Re. An extrapolation of stellar 
population models is required to account for line-strength values 
in the central region. In all panels, dashed vertical lines mark the 
seeing radius. 


sensitive features of ISpiniello et all ll2015l) . as (1) CaHl has 
been shown to be not ma tched by extended MIL ES models 
(regardless of the IMF, see lLa Barbera et aLl201^ . while (2) 
CaH2 is strongly affected by the O 2 (B-band) atmospheric 
absorption at A ~ 6900 A in our data. Concerning Ca3, the 
(observed) index tends to decrease with radius, while our 
best-fitting solutions favour a flat radial behaviour for Ca3. 
This discrepancy might be related to some radial variation 
of the [Ca/Fe] abundance ratio which is not well captured 
by our fitting approach. In fact, although we have included 
[Ca/Fe] as a free fitting parameter (see Appendix A), the 
[Ca/Fe] is mainly constrained, in this work, from the CaHK 
line, which is strongly affe cted by the abundanc e pattern of 
other alpha elements (e.g. IVazdekis et al.l 12013 ). Neverthe¬ 
less, a decreasing trend of Ca3 with radius is the opposite 
behaviour to what one would expect for a less bottom-heavy 
IMF at larger radii (as Ca3 decreases with Fb). This is con¬ 
sistent with our conclusion that the IMF of NGC1277 shows 
only a mild radial variation. In fact, as shown in Appendix A, 


CaHl index is mild (0.002 mag) further supports a rather flat 
IMF profile. 
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Figure 5. Radial IMF-slope gradient for the relic massive galaxy 
NGC1277. The left vertical axis indicates the slope of the IMF 
in the bimodal case (Fb), i.e., the slope of the IMF for stars with 
masses larger than 0.5Mq. This Fb parameter can be transformed 
into a dwarf-to-giant ratio (right vertical axi s), which is what 
the l ine-strength indices are actually tracing dLa Barbera et alj 
|20 i 3). For reference, the Milky-Way IMF value is shown as a hor¬ 
izontal dashed line. The inferred IMF gradient is mildly decreas¬ 
ing and it is very bottom-heavy up to 1.5Re. Uncertainties (the 
orange shaded region) are obtained by removing line-strength in¬ 
dices from the fitting process which might be potentially affected 
by flux calibration, telluric correction or modelling uncertainties. 
Dashed vertical line indicates the extent of the seeing disc. 



removing this index altogether from the fits does not signif¬ 
icantly change the radial trend of Fb. Regarding aTiO, the 
only point that deviates significantly from the best-fits is 
the innermost one. This could be related to some issue in 
the extrapolation of the models at high metallicity, and/or 
some peculiar abundance ratio for R=0Re. However, as the 
deviation affects only one single point, it does not change 
our general conclusions about the radial behaviour of Fb for 
NGC1277. 

The robustness of the IMF profile, against all issues 
mentioned so far, is demonstrated in Appendix A, where 
we show how the mildly decreasing behaviour of Fb persists 
under different telluric correction procedures (see §[lll , as 
well as different combinations of indices included in the fits, 
and changes in the modelling assumptions. 

Ultimately, gravity-sensitive indices trace the dwarf-to- 
giant star ratio i n the IMF, regardless o f the adopted IMF 
parametrization (iLa Barbera et aDl2013l) . This is explicitly 
shown in the right vertical axis of Fig.O where we represent 
the fraction (in mass) of stars below 0.5M© for each IMF 
slope. The inferred dwarf-to-giant ratio in NGC 1277 is, at 
all radii, at least a factor of 2 larger than in the Milky Way 
(horizontal dashed line). 

The inferred stellar population properties of NGG 1277 
(ages, metallicities, and IMF slopes; as shown in Figs. [4] 
and (SJ can be used to estimate the stellar M/L radial pro¬ 
file. Notice, however, that for old stellar systems (> 10 
Gyr), stellar population analysis cannot provide fully reli¬ 
able mass-to-light ratios, as it relie s on the assumption of a 
given parametrisation for the IMF llFerreras et al.lboi^ ). In 
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Figure 6. Radial M/L ratio (V-band) gradient for the massive 
compact galaxy NGC 1277, under the assumption of a bimodal 
IMF. At each radius, the best-fitting age, metallicity and IMF 
slope shown in Figs. [4] and [5] translate into a V-band M/L ra¬ 
tio. The right-hand axis represents the “mismatch parameter”, 
o, defined as the ratio between the inferred M/L and that ex¬ 
pected from a stellar population with the same age and metallic¬ 
ity as observed, but with a Milky Way-like IMF slope. Notice that 
the absolute values of the M/L gradie nts heavily depend o n the 
parametrisation adopted for the IMF llFerreras et al.|[2014h . The 
light shaded region accounts for the uncertainties in all the stel¬ 
lar population parameters (e.g., error bars in Fig. [4] and shaded 
region in Fig. Ell. Dashed vertical lines mark the seeing radius. 


this sense, our favoured bimodal distribu tion seems to bet¬ 
ter match the dynamical M /L estimates l|La Barbera et ^ 
I2OI3I : Fpiniello et al.ll20i3l . while a single power-law distri¬ 
bution provides unrealistically high M/L’s. For a bimodal 
IMF parametrisation, the expected M/L radial profile of 
NGC 1277 is shown in Fig. [6] 
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Figure 7. Radial Vrms profile of NGC 1277. This quantity, ac¬ 
counting for the total kinetical energy, is more closely related to 
the effective velocity dispersion, which has been found to tightly 
correlate with the observed IMF variations. The radial Vrms vari¬ 
ation is significant, albeit milder, than the cr variation. Dashed 
vertical lines mark the seeing radius. 


This value is incompatible with our inferred IMF gradient 
(AFb ~ —0.5). Note, however, that our conclusions apply 
to the local velocity dispersion, not to the effective central 
velocity dispersion. To investigate the relation between IMF 
slope and local kinematical properties of NGC 1277, we show 
in Fig. [7] its Vrms radial profile, defined as: 



4 DISCUSSION 

4.1 The excess of low-mass stars in ETGs 


What drives the excess of low-mass stars in ETGs remains 
unknown, however there are presently two distinct candi¬ 
dates, namely [Mg/Fe] and a. The latter may be viewed as 
a dynamical controlling parameter, whereas the former is a 

stellar population param eter. __ 

_ Dynamical studies dTreu et al.l 120101 : ICappellari et ahl 

l2012ll suggest that the galaxy central velocity dispersion 
correlates strongly with the IMF slope - with more mas¬ 
sive galaxies having steeper IMF slopes than less mas¬ 
sive galaxies. This result has found supp ort from sev¬ 
eral groups using different methodologies (|Ferreras et al.l 
l201,ll : iLa Barbera et all l201,ll : ISpiniello et alJ 2014h . How¬ 
ever, based on our analysis, we reject the hypothesis that 
local velocity dispersion is the principal agent driving IMF 
variations. Given the velocity dispersion gradient observed 
in NGC 1277 (from ~ 420 km s“^to 150 km s“^), if the 
correlation between IMF slope and central velocity disper¬ 
sion holds also for local velocity dispersion, we would expect 
a typical radial change in the IMF slope of AFb = —2.4. 


Note that the e ff ective velocity dispersion ana ly sed i n 
iLa Barbera et all (1201.11 1 and in ISpiniello et al.l ll2ni4h . 
which was estimated using a 3 arcsec wide fibre, is better 
represented by Vrms rather than by the local sigma, since 
the latter does not include the line broadening associated 
with the rotational velocity profile. Both Fb and Vrms de¬ 
crease with radius; Vrms is greater than 300km s“^, and Fb 
is steeper than the Milky Way value at all radii in NGC 1277. 
This is consistent with previous works reporting a bottom- 
heavy IMF at high effective velocity dispersion. 

As an alternative explanation for the IMF variations, 
based on stellar population analysis. IConrov fc van DokkumI 
(l2012tl found that the IMF slope also correlates with the 
overabundance of a-elements, in particular with the [Mg/Fe] 
ratio, in the sense that galaxies more enhanced in [Mg/Fe] 
have steeper IMF slopes. As we mentioned in ^ our [Mg/Fe] 
measurements are heavily dependent on the interpolation 
scheme we use, and we find that our data can entertain 
either slightly decreasing or increasing [Mg/Fe] radial gra¬ 
dients. Therefore, from the present study we cannot as- 
sess the relation betwe en IMF slope and ]Mg/Fe] (but see 
iLa Barbera et ah |[2^. 
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To understand the d ifferen ces among the afore¬ 
mentioned studies, ISmi t3 |I2014| has independen t ly re - 
viewed the dynamical results from ICappellari et al.1 (|2ni2l l 
and those based on stell ar population analysis of 
IConrov fc van DokkumI il2012li . Although he found an over¬ 
all agreement between both studies, the galaxy-by-galaxy 
comparison shows that the correlation between IMF slope 
and velocity dispersion does not agree with that found when 
comparing the IMF and [Mg/Fe]. 

The rapidly growing observational evidence of a vari¬ 
able IMF slope may be in tension with other well-established 
properties of ETGs in the local Universe. In particu¬ 
lar, chemical evolution models assuming a bottom-heavy 
IMF are unable to reproduce the observed metal licities 
and abundance ratios of ETGs llArrigoni et al.l [20101 1. This 
idea has been recent ly rev isited in a ser i es of articles by 
Weidner et al.l ll201^ and Ferreras et al.l (l2015l l (see also 
Vazdekis et al.l[l997l : iLarsonlllQQsh who found that a time- 


dependent IMF is necessary to explain the stellar popula¬ 
tion proper ties of massive ETGs at 2 ~ 0. Note that, ac¬ 
cording to iFerreras et al.l il2015ll . a constant-in-time IMF 
can be ruled out, not only for a bimodal IMF shape, but 
al so for a two-segment IMF para metrization (e.g. as defined 
in iGonrov fc van DokkumI |2012| 1. and even for a V-shaped 
IMF (characterized by an enhanced fraction of both low- 
and high-mass stars). 


4.2 The emerging picture 

Our analysis of the stellar population properties of 
NGC 1277 is summarized in Figs. |T]and[S] The radial age 
profile of this galaxy, even allowing for possible IMF varia¬ 
tions, indicates that the underlying stellar population is very 
old (> 10 Gyr) at all radii, w ith no evidence of recent star 
formation (iTruiillo et al.l[2014h . This result suggests that the 
stellar component of NGC 1277 was formed at high redshift 
{z > 2), and that solving for the IMF at each radius does not 
convert NGC 1277 into a young, massive and compact sys¬ 
tem, as those identified in the local Universe (ITruiillo et al.l 
I2OO9I: iFerre-Mateu et al.ll2012l, |2013|) . 

In terms of kinematics, NGC 1277 shows an extreme 
radial velocity dispersion gradient, ranging from more than 
400 km s“^in the centre to cr 150 km s“^beyond 1 Re. 
The rotation in this galaxy shows the opposite behaviour, 
rising to V)ot ~ 280 km s“^at 1 Re. Interestingly, these dif¬ 
ferences in the kinematical properties between the centre 
and the outskirts of NGC 1277 are not reflected in the 
properties of its stellar populations. The radial variation 
of a-element overabundance, age and IMF slope are mild, 
whereas the total metallicity decreases by ~ 0.4 dex from the 
centre to 1 Re. These gradients ar e consistent with a fast, 
mono Uthic-like formation process lISanchez-Blazauez et al.l 
[20 o 3). In this sense, in terms of its stellar population prop¬ 
erties, NGC 1277 does not differ from typical massive ETGs 
in the Local Universe. However, its extreme dynamics and 
compact morphology set it apart from the vast majority of 
local ETGs. 

The IMF radial profile of NGC 1277 further constrains 
the formation process of massive ETGs. NGC 1277 shows 
much weaker radial variations of the IMF slope (AFb = 
—0.5) when compared to th e IMF gradient of t h e two 
massive galaxies analysed by iMartfn-Navarro et al.l (I 2 OI 5 I I 


(AFb = —1.5 ). Our suggestion is that NGC 1277 repre¬ 
sents the monolithic-like formation phase of massive ETGs. 
In this picture the later accretion of less massive satellites 
(with “standard” Kroupa-like IMF slopes), which would set¬ 
tle preferentially in the galaxy periphery due to their lower 
densities, would steepen the pristine IMF gradient mak¬ 
ing i^jnore_similm_to_^atjound in nearby massive galax¬ 
ies [MartfmN^^rro^^r^^ ll2015h . 

Indeed, we hypothesize that relic galaxies, such as 
NGC 1277, will evolve to be come the cores of m assive ellipti¬ 
cals in the nearby Universe (ITruiillo et al.l2014ll . The bulk of 
the stellar mass of this galaxy was formed in a short (r ~ 100 
Myr) and intense (SFR > 10® MQ/yr) burst at redshift 2 > 
2, but remained concentrated within a few kpcs in the centre 


of the dark matter halo. Afterwards, through dry merging 

( van Dokkum et ^l2010l : ITruiillo et al.llioilt IFerreras et al.1 

l2014ll . the system would typically evolve to the charac¬ 
teristic sizes and stellar masses that are observed at 2 ~ 
0. This in terpretation is supported by numerical simula - 
tions (e.e. lOser et al]l201ol : iNavarro-Gonzalez et al.ll2013l L 
and also by s tellar population studies of the outskirts of 
nearby ETGs (ICoccato et al.l 2010j^ La^arbera et aDl2012l : 


IPastorello et al.ll2014l : iMontes et al. l2014l L Since dry merg¬ 

ers cannot significantly alter the properties of the stellar 
populations, a relic massive galaxy should present the same 
stellar population properties as the core of a typical 2 ~ 0, 
massive ETG, but with the latter’s morphology and kine¬ 
matics shaped by successive mergers. 

It may seem surprising, at face value, that such a relic 
system should be identified in such a dense and dynamic 
environment such as that re presented by the Pe rseus clus¬ 
ter (M 200 = 7.74 X 10^'^Mm. lAleksic et al.l[201oh . However, 
this is in fact a natural expectation in the hierarchical clus¬ 
tering paradigm since overdensities in what are now galaxy 
clusters were th e first structures t o collapse and to form 
massive galaxies (IMo fc Whitel[l99^ . In addition, the active 
intracluster environment may have inhibited the subsequent 
formation of younger stellar populations within NGC 1277 
either via the stripping or heating of cold gas (|Edg eet__alJ 


eitner via tue stripping or neatmg 01 cola gas (mqg eet__aij 
ll992l : [^Dea et aIll200a L In this sense. Istringer et al.l ( 201 
have recently shown tha t in the BOLSHOI numerical simula¬ 
tion dKlypin et al.llioilh . dark matter haloes similar to those 
hosting NGC 1277-like objects have experienced below- 
average mass accretion, or even mass loss, since 2 = 2. 

Finally, we find that the radial IMF behaviour seen in 
NGC 1277 is not driven by the local velocity dispersion, 
and therefore, it cannot explain the IMF variations seen 
in ETGs. The fact that the IMF slope c an vary radially 
within galaxies (iMarti'n-Navarro et al.l[2015f) suggests a com¬ 
plex galaxy formation scenario, with significantly different 
enrichment and supernovae-feedback rates when comparing 
the centre and the outskirts of ETGs. Thus, more observa¬ 
tional efforts are needed to better understand the IMF be¬ 
haviour in nearby galaxies, investigating what stellar popu¬ 
lation properties can explain the observed IMF variations. In 
this s ense, large resolved s pectroscopic surveys such as CAL- 
IFA (ISanchez et al.|[2012l l and MaNGA will provide statis¬ 
tically significant samples. In addition, the joint analysis of 
the IMF inferred from dynamical modelling and from stellar 
population analysis can reveal important information about 
the interplay between galaxy stellar population properties 
and dynamics. The intrinsic degeneracies of these methods 
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{M/L ratios, dark matter content, stellar remnants and cen¬ 
tral black-hole mass) may be broken if both approaches are 
combined. 
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APPENDIX A: ROBUSTNESS OF THE 
INFERRED IMF PROFILE 

To check the reliability of the IMF gradient of NGC 1277 
(Fig- Ell, we have explored three possible sources of systematics, 
namely, errors on telluric absorption correction, the effect of non¬ 
solar elemental abundances, and the nebular emission contamina¬ 
tion to the line. These tests — supporting the robustness of 
our results - have been performed because of some discrepancies 
between the best-fitting and observed line-strengths (see Sec. o, 
in particular regarding the Ca and Na spectral lines (see 

panels (c), (i) and (j) of Fig. EJ. 

The fact that the best-fitting lies above the observed 
values (panel (c) in Fig. 0 may indicate some residual contam¬ 
ination from nebular emission to the line profile. In the present 
work, the nebular emission correction to has been estimated 
in a self-consistent manner (see La Barbera et al. 2013 for details), 
i.e., for each IMF slope (Fb) explored in the fitting process, we ap¬ 
ply the correction calculated by simultaneously fitting the line 


profile with a 2 SSP model (for a given Fb) plus a gaussian emis¬ 
sion line. For the two innermost radial bins, the correction turned 
out to be fully consistent with zero regardless of Fb, while in the 
bins from 0.5 to 1.5 Re, a mild correction was foundQ, in the 
range from 0.1 to 0.2 A. Therefore, it is unlikely that differences 
in panel (c) of Fig. 0 are due to residual nebular contamination. 
Note also that underestimating (in particular for the outer¬ 
most bins) would imply that, given the age sensitivity of the TiO 
features, the inferred IMF slope is shallower than the “true” one, 
hence strengthening our conclusion about a shallow IMF gradient 
in NGC 1277. The possible explanation for the mismatch between 
the best-fitting and observed may reside in the uncertainty of 
the stellar population models used to account for non-solar abun¬ 
dance ratios on line strengths (in particular H^^). We modelled 
deviations f rom solar partition with the semi-empirical approach 
described inlLa Barbera et al.l ll2013h . with the further aid of the 
IConrov & van ^okkuml l|2012l . hereafter CvD) stellar population 
models. In practice , the terms in Eg. [3 ca n be expanded as 

follows (see §3.3 in lMartin-Navarro et al1]2015lb 

■ [a/Fe] -h Ax,i ■ SX 

where Ca,i is the observed sensitivity of the i- th index to [a/Fe] 
(deri ved from SDSS stacked spectra, as in iLa Barbera et ^ 
l2013fl : Ax, i is the theoretical sensitivity of each index to the 
abundance ratio of a specific element X, estimated with CvD 
models, and 5X is a “residual” abundance of element X, not ac¬ 
counted for by the Coc,i correction, and treated as a fitting pa¬ 
rameter. In other terms, we remove the bulk dependence of each 
index on [a/Fe] with the Ca^i terms, while still allowing for some 
residual abundance effect o n the indices through t he best-fitting 
5X’s. In practice, following iLa Barbera et al.ll201^ . we have con¬ 
sidered only the effect of elemental abundances of selected gravity- 
sensitive features, i.e. X = Ca, Na, and Ti. Regarding we 

have AX ~ 0 for all X’s but X = Ti, i.e. CvD models predict 
a significant sensitivity of to [Ti/Fe] abundance ratio (in 
the sense of increasing with [Ti/Fe]). In fact, we found that 
the [Ti/Fe] tends, on average, to increase the best-fitting by 
~0.15A, possibly explaining at least part of the observed mis¬ 
match between the best-fitting 0 and observed in panel (c) 
of Fig. [3] Since the dependence of on specific elemental abun¬ 
dances may be significantly affe cted by theoretical uncertai nties 
on stellar population models ('see lCervantes fc Vazdekij|2009l) . we 
have repeated the fitting process by (i) removing the ATi,H^^ -^Ti 
term from Eq.[3] and (ii) excluding the index altogether from 
the fits. The resulting IMF profiles from these tests are plotted 
in Fig. rm showing no significant differences with respect to the 
reference profile in Fig. [5] Moreover, neglecting the [Ti/Fe] ef¬ 
fect on leads to a better agreement between observed and 
best-fitting radial trends for this line, as shown in Fig. \A2\ 

As mentioned in § [2] one major issue in constraining the 
IMF is the removal of telluric lines from the observed spectra. 
In the present work, we have removed atmospheric absorption 
with two independent approaches, i.e., (1) using a spectropho- 
tometric standard star and (2) applying the software MOLEC- 
FIT. Fig. \M\ compares the corresponding IMF best-fitting pro¬ 
files, showing no significant differences within the uncertainties. 
The two approaches lead to an estimated IMF radial variation of 
(1) AFb = —0.5 and (2) AF^ = —0.6, respectively. 

To further explore to what extent the mismatch of Ca and 
Na lines for some radial bins (see panels (i) and (j) of Fig. 3) 
may be affecting our conclusions, we have recomputed the IMF 


Notice that this correction is already applied to observed 
values (blue dots) in Fig. [3] 

^ Note that the best-fitting indices in Fig. [3] are actually given 
by EWM,i + Ac,,i. 


The IMF of a massive relic galaxy 9 



S 

CO 

O) 

O) 

aj 

E 

§ 

o 


o 

c 

o 


o 

CO 


Figure Al. Same as Fig. [5]but including tests to check the ro¬ 
bustness of our IMF inference. The gray solid line indicates the 
IMF profile obtained using Molecfit to correct telluric absorption, 
whereas the gray dashed line corresponds to the best-fitting so¬ 
lution if the [Ti/Fe] dependence of ignored. In solid black, 
doted and dashed lines we show the IMF profile when excluding 
from the fitting the Calcium, Sodium and features, respec¬ 
tively. The overall agreement among all the fits shows that our 
main result is not significantly affected by telluric, nebular emis¬ 
sion or model uncertainties. Dashed vertical lines mark the seeing 
radius. 



Figure A2. Observed and best-fitting radial profiles. The 

right panel is the same as panel (c) in Fig. (3] whereas the left 
panel corresponds to the case where [Ti/Fe] abundance effects on 
are neglected. Notice, in the latter case, the better agreement 
between best-fitting and observed indices, with no significant im¬ 
pact on the inferred IMF profile (Fig. lAlb . 


profile by excluding, as we did for (see above), all Na and Ca 
lines in turn, from the fitting procedure. Again, the resulting IMF 
gradients, shown in Fig. [ah are compatible with the IMF profile 
reported in Fig. [5] It is worth noting that the impact of removing 
Na lines is the most significant, suggesting either some residual 
telluric contamination in the data, or some model uncert ainties 
in dealing with the abundance ratios ijSDiniello et al.ll2015h . 

All the above tests, summarized in Fig. rm indicate that our 
IMF inference is robust against i) nebular emission correction, w/, 
abundance effects (on line). Hi) age determination through 
itself, iv) telluric correction, and v) the set of lines included 
in the analysis. 











































